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I. INTRODUCTION 



Since the discovery of quasi-solvability in one- dimensional quantum mechanics [l| , many 
new ideas and concepts have been discovered and developed. The progress until the early 
1990s is summarized in Ref. So far, most of the known quasi-solvable quantum systems 
of one-degree of freedom are related to the ones constructed from the 5l{2) generators j^. 
This family of quasi-solvable models was completely classified under the consideration of 
the GL(2,M) invariance of the models Recently, a new type of orthogonal polynomial 

systems was discovered in connection with a quasi-solvable model 0]. Soon later, it was 
shown that similar type of polynomial systems, known as weakly orthogonal polynomial 

In the last 



systems, can be constructed for all the sl{2) quasi-exactly solvable models 
couple of years, much attention has beenpaid to quasi-solvable models from the viewpoint 
of TV-fold supersymmetry B [11 [ll|, [11 El H H H H H H H Furthermore, 
several reahstic physical systems have been found, which can be reduced to one-dimensional 
quasi-solvable models HHIlllilSHIiiliiHIsllli. 

Unlike one-body systems, construction of quasi-solvable quantum many-body systems 
suffers from another difficult problem. The difficulty originates from the fact that second- 
order differential operators of several variables are, in general, no longer equivalent to the 
Schrodinger operators. Indeed, construction of quasi-solvable two-body problems by the 
naive extension of the 31(2) method to the rank two algebras, investigated extensively in 



the early 1990s 0, M, IH, 0, 37, [sll, led to the Schrodinger operators in curved 

space. In 1995, a breakthrough was achieved in Ref. [s^, where the exact solvability of the 



rational and trigonometric A type Calogero-Sutherland (CS) models |40l . 41, i^l for any 
finite number of particles were shown by a similar algebraic method. The key ingredient is 
to construct the generators of sl{M + 1) Lie algebra in terms of the elementary symmetric 
polynomials which reflect the permutation symmetry of the CS models. This idea was 
further employed to show the exact solvability of the rational and trigonometric A and BC 
type CS models and their supersymmetric gene ralizations (43l |. and to show the quasi-exact 
solvability of various deformed CS models ji^. l45l |. So, it is natural to ask what kinds of 
quasi-solvable quantum many-body systems can be constructed from the same algebraic 
procedure. Recently, this classiflcation problem was partly accessed in Ref. )46| though the 
method depends on the speciflc ansatz employed. In our previous short letter |47|, we have 
briefly reported a partial answer for this problem without recourse to any speciflc ansatz 
by restricting the consideration up to two-body interactions. In this paper, we will give the 
details and the further developments on the issues. Especially, we will see that the GL{2, K) 
symmetry plays a central role not only on the classiflcation problem but also on the problem 
of what kind of many-body operators can exist in the quasi-solvable quantum Hamiltonian. 

The paper is organized as follows. In the next section, we give the deflnition of quasi- 
solvability together with its important subclasses such as quasi-exact solvability, solvability 
and so on. Based on the deflnition, we present a procedure to construct a quasi-solvable 
model by introducing a set of differential operators. In Section IIIIl we show that, by the 
introduction of elementary symmetric polynomials and suitable choice of the set of differ- 
ential operators which determines a solvable sector of the model, a family of quasi-solvable 
models can be constructed from a polynomial of the sl{M + 1) generators of second- degree. 
Section IIVI deals with the canonical-form condition, namely, under what conditions we can 
obtain the quasi-solvable quantum systems from the algebraic operators constructed in Sec- 
tion mil We explicitly show the procedure to solve the canonical-form condition for the 
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two-body problem. In Section we explain the GL{2, K) symmetry of the model. It turns 
out that this symmetry is so restrictive that the quasi- solvable quantum two-body model 
obtained in Section HVl can be uniquely generalized to M-body (M > 3) systems. An inter- 
esting relation between (quasi-) solvability and supersymmetry is discussed. In Section IVH 
we show that the many-body quantum systems obtained solely by the consideration of the 
GL{2, K) symmetry are actually quasi-solvable. This is achieved by expressing the models 
in terms of the elementary symmetric polynomials. In Section IVIH we classify the mod- 
els under the consideration of the GL{2, K) symmetry, by essentially the same way as in 
Ref. 01 ■ The explicit forms of both the superpotentials and potentials are presented. In 
Section rVIIIl we investigate the possibility of including M-body (M > 3) interaction terms 
without destroying the quasi-solvability. Again, the GL{2, K) symmetry plays an essential 
role in this problem. We give several concluding remarks in the final section. Some useful 
formulae are summarized in Appendices. 



II. QUASI-SOLVABILITY IN MANY-BODY SYSTEMS 

First of all, we will giv e the definition of quasi-solvability and some notions of its special 
cases based on Refs. (igL l38| . A linear differential operator H of several variables q = 
(gi, . . . , qm) is said to be quasi-solvable if it preserves a finite dimensional functional space 
Vj\f whose basis admits an explicit analytic form: 

HVj^cVjv, diraVj^ = n{J\f) < oo, Vat = span {0i(g), 0„(^)(g)} . (2.1) 

An immediate consequence of the above definition of quasi-solvability is that, since we can 
calculate finite dimensional matrix elements Sk,i defined by, 

n{Ar) 

H(t>k=Y,Sk,i(t)i (A; = l,...,n(Ar)), (2.2) 

1=1 

we can diagonalize the operator H and obtain its spectra on the space Vv, at least, alge- 
braically. Furthermore, if the space VV is a subspace of a Hilbert space L'^{S) {S C M*^) 
on which the operator H is naturally defined, the solvable spectra and the corresponding 
vectors of Vj\f give the exact eigenvalues and eigenf unctions of H, respectively. In this case, 
the operator H is said to be quasi- exactly solvable (on S). Otherwise, the solvable spectra 
and the corresponding vectors of Vj^f only give local solutions of the characteristic equation 
and have, at most, restrictive meanings in the perturbation theory defined on the physical 
space L\S) BHII^. 

A quasi-solvable operator H of several variables is said to be solvable if it preserves an 
infinite fiag of finite dimensional functional spaces Vv, 

Vi C V2 C ■ ■ ■ C Vat C ■ ■ ■ , (2.3) 

whose bases admit explicit analytic forms, that is, 

HVj^ C Vat, dim Vat = n{Af) < 00, Vx = span {0i(g), . . . , (j)n{Af)iQ)} , (2.4) 

for Af = 1, 2, 3, . . .. Furthermore, if the sequence of the spaces (j2.3p defined on 5 C M-^ 
satisfies. 



VuiS)^L\S) (AT ^00), (2.5) 
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the operator H is said to be exactly solvable (on S). 

In the case of a system of a single variable, it was proved that quasi- solvability is essentially 
equivalent to Af-fold supersymmetry In this context, it is convenient to formulate the 
quasi-solvability in terms of A/'th-order linear differential operators. Let us introduce A/'-fold 
supercharges by, 

Qa/- = E * = E (2.6) 

{i} {•«} 

where ip^ij and V'jj} are fermionic coordinates and Pj}^ are AAth-order linear differential 
operators of several variables g^. Using these operators, we define the vector space V^f as, 

V^ = f|kerP«. (2.7) 

Now, it is easy to see that an operator H is quasi-solvable with the solvable sector (j2.7|) if 
the following quasi-solvabihty condition holds 

P^^HVm = VW. (2.8) 

Therefore, the above formulation gives a concrete procedure to construct a quasi-solvable 
operator. In practice, however, the problem often becomes more tractable when we make 
a suitable similarity transformation on the vector space and on the operators and make a 
suitable change of variables. In the case of M = 1, the procedure was successfully employed 
to construct the type A A/'-fold supersymmetric models ISl- following section, 

we will generalize the procedure to the case of several variables. 



III. sl{M + 1) ALGEBRAIZATION OF THE MODELS 

Let us consider an M-body quantum Hamiltonian for a system of identical particles on 
a line, 

I ^ g2 

Hm = -^^-g^ + y{(li,---,(lM), (3.1) 

i=l * 

which possesses permutation symmetry, that is, 

V{...,qi,...,qj,...) = V{...,qj,...,qi,...) "rfi^j. (3.2) 

To algebraize the Hamiltonian (j3.H) . we will proceed the following three steps. First, we 
make a gauge transformation on the Hamiltonian (|3.1|) : 

HM = e^^''^H^e-^^''l (3.3) 

The function VV(g) is to be determined later and plays the role of the superpotential when 
the system Eq. ()3.1|) is supersymmetric. As in Eq. ()3.3p . we will hereafter attach tildes 
to both operators and vector spaces to indicate that they are quantities gauge-transformed 
from the original ones. In the next step, we change the variables qi to hi by a function 
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/i of a single variable; hi = h{qi). Note that the way of changing the variables preserves 
the permutation symmetry. The third step is the introduction of elementary symmetric 
polynomials of hi defined by, 



M 



au{h)= J2 h^,■■■h,, {k=l,...,M), (3.4) 



ti<---<ik 



from which we further change the variables to a^. Owing to the permutation symmetry, 
the gauged Hamiltonian ()3.3|) can be completely expressed in terms of these elementary 
symmetric polynomials ()3.4|) . We choose a set of components of the AA-fold supercharges in 
terms of the above variables ai as follows: 

^i^^ = 7, ^ (l<^i,...,W<M), (3.5) 

dcTi, ■ • ■ daij^ 

where {i} is an abbreviation of the set {ii, . . . ,ij\f}. Using these A/'-fold supercharges, we 
define the vector space by Eq. ()2.7|) . in which, 

Vjsf = span 1^;^^ ■ ■ ■ alf : rn G Z>o, 0<^ni<7V-l|. (3.6) 

For given M and A/", the dimension of the vector space ()3.6|) becomes, 

J. -^'(n + M-l)! {U + M-l)\ 



n=0 



We will construct the system ()3.3|) to be quasi-solvable so that the solvable subspace is given 
by just Eq. ()3.6p . This can be achieved by imposing the quasi-solvability condition of the 
gauge-transformed form [3, 12 1, 

P«^ArVAr = V{z}. (3.8) 

The above condition ensures Hj^Vj^ C Vat, that is, the subspace is invariant under 
the action of Hj^. As a consequence, we can diagonalize the gauged Hamiltonian ()3.3|) 
in the subspace fl3.6|l with the finite dimension given by Eq. ()3.7|1 . and can obtain the 
corresponding spectra algebraically. Then, it is evident that they are also the spectra of the 
original Hamiltonian (j3.1|) in the subspace Vv given by, 

Vm = span i.a'l'--- cx^f e^^^^^ : Hi G Z>o, 0<^r2, <Ar-l|. (3.9) 



i=l 



Thus, the system Hj^ turns to be quasi- solvable. 

The general solution of Eq. ()3.8p can be obtained in the same way as in Ref. 17[. As in 
the case of the one-body models, it is sufficient to find differential operators up to second- 
order as solutions for Hj\f since we are constructing a Schrodinger operator in the original 
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variables g^. We can find that solutions expressed in terms of first-order derivatives are the 
following: 

A.Eo. (3.10a) 
= Ei,, (3.10b) 



These + 2M operators -E'ka(/« + A > 0) together with Eqq constitute the gl{M + 1) Lie 
algebra: 

[Ekx, Ef^iy] = S^^xEku — 5k.,uE^x- (3.11) 

Here we note that throughout this paper the Latin indices take values from 1 to M while 
Greek indices from to M, as in the standard convention in relativistic theories. Solutions 
involving second-order derivatives are of the following form: 

E^,Eo,, (3.12a) 



daidaj 



dakdai 

M 

a. ' 



Q2 

^i-Q^-g^^ = Eij^ok, (3.12b) 
— a — ^ EikEji — 5j^kEii, (3.12c) 



^ - 2 - E ^''d^ - 1 - E 5^ = ^^o^^o- (3.12e) 

k=l I'/ \ i^i U 



Therefore, the general solution which contains up to the second-order differential operators 
reads, 

M M 

Hj\f = — ^K\,^luE^,\E^^ + ^ B^xEkx — C, (3.13) 

where Af^x,fiu, Bi^x, C are arbitrary constants. Since Eqq + X]f=i En = Af — I, the algebra 
represented by Eqs. (j3.10p is essentially s[(M -|- 1) and we can set ^^^,00 = Aoo,fj.u = 
Vk, A,/i, u and Bqq = 0. However, all the terms in Eq. (|3.13|) are not linearly independent 
and thus we can set some of the constants zero. From the Lie algebraic relations (j3.1ip . we 
can put, 

(3.14a) 
(3.14b) 
(3.14c) 
(3.14d) 



AiOJO — 


■ ^iOJl 


= 


Vi < J, 




^ikjl 


= 


Vi < j or VA; < / 


Aikfil = 




= 


< I, 


Aok,jl = ^OfcJO = 


Aik,jO 


= 
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Furthermore, we can remove the remaining redundant terms as follows: 

M 

EioEoi = -Y,E^kEki + {M-Af+l)Ea A,o,oz = Vi,/. (3.15) 



k=l 



After the substitutions (|3.14p and (|3.15p . the remaining terms in Eq. (j3.13|) turn to be linearly 
independent. Substituting the generators (j3.1Up for Eq. (|3.13|) and taking Eqs. (j3.14p and 
(|3.15p into account, we obtain the general solution (|3.13p in terms of af 



Q2 



k,l=l 
M 



dakdai 



+ J2 [Bo{(r)ak - Bk{a)] - C{a), (3.16) 

k=l ^'^ 

where A^, Aki, -Bk and C are the following second-degree polynomials of several variables: 



M 

^k(c^) = y^^Aoj^cTjCrj, (3.17a) 

M M 

Aki{(y) = ^ Aikjiaiaj + ^ Aik^oiai + Aokfii ik>l), (3.17b) 

i>j i=l 

AI 

Bo{a) = 2{Af - 2)Ao{a) - BioCTi, (3.17c) 

i=l 

M M 

Bk{a) = {X- 2)Ak{a) + ^ A^jk^i - ^ Bikai - Bok, (3.17d) 

i>ji>k) i=i 

M 

C{a) = {^/- l)(Ar-2)Ao((T) - (AT - 1) ^ 5,o(Ti + a (3.17e) 

i=l 

Before closing the section, we will investigate the condition under which the quasi-solvable 
operator ()3.16| becomes solvable. If we act the sl{M + 1) generators for a fixed M on the 
vector spaces Vk for any A; = 1, 2, 3, . . ., we have, 

EoiVkCVk-i, EijVkCVk, EMcVk+i. (3.18) 

Thus, the operator ()3.13|) preserves an infinite flag of the spaces Vk given by Eq. ()3.6|) . 

ViCV2C---CVfeC--- , (3.19) 

if it does not contain the operators such as EioEjQ, Ei^Ejk and EiQ at all. If we note from 
Eqs. fl3l4g|l and (ITTH|) that we can arrange the operator (j3.13p so that it does not contain 
operators like Eq^E^q and EiqEqi, we conclude that the operator is solvable if it does not 
contain the raising operator Eio at all. Thus, the solvability condition can be written as, 

Aioj^ = B,o = Vz,j,/€. (3.20) 
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Finally, if we substitute the above condition for the set of polynomials ()3.17|) . we get the 
solvability condition for the gauged Hamiltonian ()3.16|) as, 

A^{a) = Bo{(j) =0 Vk. (3.21) 

If this is the case, Sfc(cr)'s become polynomials of at most first-degree while C(cr) becomes 
a constant: 

M M 

Bk{<y) = ^ AijjkO'i — ^-BifeCTj — Bok, (3.22) 

i>j{>k) 1=1 

C{a) = C. (3.23) 



IV. CANONICAL-FORM CONDITION 

It is evident from the construction that any operators obtained from Eq. ()3.16|) by gauge 
transformations are quasi-solvable. However, this does not necessarily fulfil our purpose at 
the quantum level. One of the most difficult problems one would encounter in the algebraic 
approach to the quasi-solvable quantum many-body systems is to solve the canonical-form 
condition: 



= e-^(^)^^e>-(') = -IJ:^ + V{q). (4.1) 



2^ dq: 

i=l ^' 

If the Hamiltonian ()3.16p is gauge-transformed back to the original one, it in general does 
not take the canonical form of the Schrodinger operator like Eq. (13.11) and one can hardly 
solve, for arbitrary M, the conditions under which a gauge-transform of Eq. ()3.16jl could 
be cast in the Schrodinger form. This difficulty can, however, be overcome in our case by 
considering the symmetries. 

Suppose we can solve the canonical-form condition for an M. Then we would obtain a 
quasi-solvable Hamiltonian which might contain up to M-body interactions: 

i=l * 2=1 i<j 

M 

■■■ + g„,{M) Vm{qh,---,qiJ + --- + 9M{M)VM{qi,...,qM), (4.2) 

ii<---<im 

where each gm{M){l < m < M) denotes the coupling constant of the corresponding m- 
body interaction and may depend on the number of the particle M. Then, we can get a 
quasi-solvable M-body model with up to M'-body interactions (M' < M) if we turn off all 
the coupling constants g^iM) for m > M': 

Hm= -^J2d^ + 9i{M)Y + 92{M) Y V2{q^, qj) + ■■■ 

2=1 2=1 i<j 

M 

... + g^,{M) Yl yM'{q^,,...,q^,„). (4.3) 

2l<--<2,,/ 
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The resultant model should be, when we put M = M', one of the models constructed 
from the si{M' + 1). This is a consequence of the permutationally symmetric construction. 
Actually, the gauged Hamiltonian ()3.13p composed of the s[(M + 1) generators is reduced 
to the one composed of the generators of the subalgebra s[(M' + 1) C sl{M + 1) if we set 
hi = ioT i > M' involved in both the generators (j3.1(J|) and the representation space (j3.6|) . 

The above observation means that, we can know the functional form of the M-body quasi- 
solvable potential which involves up to M'-body interactions (M' < M) if we can solve the 
canonical-form condition for the M'-body case. Furthermore, as we will later see in section 
IVl the M-dependence of each coupling constant gm{M) can be uniquely determined by 
considering another symmetry, namely, GL{2, K) symmetry. Therefore, it is necessary and 
sufficient to solve the M'-body canonical-form condition as far as we are concerned with the 
quasi-solvable models containing up to M'-body interactions. The simplest but non-trivial 
case is the two-body case. In the followings, we will see that we can actually solve the 
canonical-form condition for the two-body problem. 



A. The Conditions for the Two-body Case 

To solve the condition, we must first transform the variables from cTj back to hi. For this 
purpose it is convenient to note that hi and h2 are the solutions of the algebraic equation: 

h1-aihi + a2 = Q (^ = 1,2). (4.4) 

From the differential of the above relation: 

d/z, = ^, ^' dai - — i da2 (z = 1, 2), (4.5) 

Zhi — ai 2hi — ai 

we can easily obtain, 

^ = yJ^ ^, A = _y-^ ^. (4.6) 

dai ^ hi — hi dhi ' da2 ~^ hi — hj dhi 

Substituting Eqs. (j3.4j) and (j4.6|) for the gauged Hamiltonian (I3.13|) . we can express it in 
terms of hi as follows: 



1 ^ 

- (h \2 Yl [^0^^ ~ ^1^*^ + + Au)hl - A2ihi + A22] 



\ ^o [2Ao(T^ - Aia20i + A2al - 2(^2 - An)(72 

[hi - h2y 

, I f d' 1 d 

- A2iai + 2A22J ( + ^ 



, ^ Bph^ - BiK + B2 d 



. , . -''I ""I 
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Let us first consider the canonical-form condition with respect to the second-order differential 
operators. Recalling the relations: 



dhi h'i \dqi dqi 



-w 



92 



dhidh2 



(92 



h[h'2 ydqidq2 ' ^ dqi dqj dqidq2 dqi dq2 ) ' 

row _hi\d_ r ow 

dqf V 9qi h[ J dqi dqf \ dq^ 



h'ldW 
h'i dqi 



(4.8a) 
(4.8b) 

(4.8c) 



and requiring the second-order differential operator to be gauge-transformed back to the 
Laplacian in the flat space, we obtain the following conditions: 



2Ao(J2 - Aio-20-i + ^20-? - 2(A2 - Aii)o-2 - A210-1 + 2A22 = 0, 

1 



(4.9) 



A,ht - A,h^ + (A2 + An)h', - A2ih, + A22 = -{h - h2)\K)\ (4.10) 



Once the above conditions satisfied, the gauged Hamiltonian is transformed to. 



i=l 



^ raw _hl\d_^ ^ /dW 
dqf V % h'J dqi dqf \ ^q^ 



h'ldW 
h'i dqi 



^ ^ Bphj - B,h, + B2 ^ d ^dW 



{hi - hj)h\ 



dqi dqi 



C. 



(4.11) 



The condition that the above operator does not contain the first-order differential operator 
reads. 



dW h'l , Boh]-Bih, + B2 

+ ^-rv-, (^7^J, 



(4.12) 



dqi 2h[ {hi - hj)h[ 

Finally, with the condition ()4.9|1 . ()4.10j) and ()4.12|1 satisfied, the original Hamiltonian takes 
the desirable Schrodinger form as. 



e-^^^e^ 



d' 1 ^ 

"2^5^^ 2^ 

i=i ^* i=i 



V dqi ) dqf 



C{a{h)). (4.13) 



In the following, we shall solve the obtained canonical-form conditions ()4.9|1 . ()4.10j) and 
()4.12|1 in order. 



B. The first condition 



The first condition ()4.9j) is an algebraic identity and thus is satisfied if and only if all the 
coefficients of the polynomial vanish, i.e.. 



^20,20 


= 0, 


2^20,10 


— ^20,21, 




2^10,10 


— ^20,11 "~ ^20,22; 


^10,11 


= ^20,12, 


^10,12 


= 0, 




^20,22 


= ^21,21, 


2^20,12 


= 2y42l,ll — ^22,21, 


2v4io,i2 


= 2An,n 


— ^22,11, 


^12,11 


= 0, (4. 


^21,01 


= —^22,22? 


^22,01 


= 2v4ii^oi 


+ 2v422,12! 


^12,01 


= 2^12,12, 


^01,01 


= —^22,02, 


^02,01 


= 2Ai2,02 




^02,02 


= 0. 
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There are 27 parameters A^x^^y in Eq. ()3.16p for M = 2. The above conditions reduce the 
number of these free parameters to 12. 



C. The second condition 

With the first condition ()4.14|) satisfied, the l.h.s. of the second condition ()4.10|) becomes, 
AX - Aihl + (A2 + An)h^ - A^iK + A22 

= (^20,10O-2/if + ^10,100-1 /if - A20,22O"2/if + ^21,110^2/^^ + Au^uCTlh'^i 

-A22,22CT2hi + Aii^Qihl - A22,Vl<y2 " ^22,02/ii - ^12,120-1 - ^12,02) {^i " hj) 



- A2o,i2h^ih, + 2hj){hi - hj 



(4.15) 



The differential equation ()4.10j) has a solution if and only if Eq. ()4.15|) is divided by the 
factor [hx — h2)'^. This factorization condition reads. 



^20,10 
^22,22 



0, 

2An,ii, 



^20,22 
^22,12 



2^10,10? 

^11,01; 



^21,11 
^22,02 



3A 



20,12; 
2^12,12; 



A 



12,02 



0, 



(4.16) 



from which the number of the free parameters A^^x^^^ is further reduced to 5. Then, the 
second condition turns to be, 



i\2 



AiQ^i^hl — A2Q^i2h\ + Aiiiih^ + AiiQihi + v4i2,i2 



= a^hl + a-ihl + 02/1^^ + aihi + oq = P{hi). (4-17) 

Under the conditions ()4.14|) and ()4.16p . the polynomials and A^i are expressed in terms 
of the five parameters (p = 0, . . . , 4) as follows: 

4a4(j2ai — flso"!, 

2a4(J2 ~ '^Z^2^li 

2a4(T2 — 3a3(T2cri + a20"^ — 2a20"2 + aiCi + 2ao, 



Ai(a 
A2{a 
Aii{a 
-421 (cr 
^22 (f^ 



-Aa^al + 2a2a2(Ti + 4ai(72 + 2aocri, 
2a20"2 + ciicr20"i + c^o'^i ~ 2aoO"2- 



(4.18a) 
(4.18b) 
(4.18c) 
(4.18d) 
(4.18e) 
(4.18f) 



D. The third condition 



The third condition (j4.12|) is a set of partial differential equations with two variables. In 
order that it has a solution, we must impose the integrability condition: 

d dW d dW 

(4.19) 



Using Eq. (HTT^ . we get. 



d dW 
dqj dqi 



dq2 dqi dqi dq2 



_d_ 

dqj L 



h" 



Boh^ -B,h, + B2 



{hi - hj)h'- 



(4.20) 
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where (p = 0, . . . , 4) are given by, 



Co = —Bq2, ai = Bqi — 2Bi2, 

02 = 2Bn - B22, % = i2Af + 3)03 - 25io + B21, 

04 = 2(7V - 2)04 - 520. 

Therefore the integrabihty condition (j4.19p is equivalent to, 



(4.21) 



(4.22) 



It is evident that the above relation holds if and only if both of the quantities are a constant, 
—2c. Finally, the integrabihty condition reads. 



(p = 0,...,4). 



(4.23) 



Combining Eqs. ()4.21|1 and ()4.23|1 . we obtain. 



B02 = 2cao, Boi - 2Bi2 = -2cai, 

2Bu - B22 = -2ca2, 25io - B21 = {2Af + 3 + 2c)a3, 

52o = 2(Ar-2 + c)a4, 



(4.24) 



from which the number of the free parameters B^x reduces to 3. When the integrabihty 
condition is fulfilled, the third condition (j4.12j) turns to be. 



dW h'! 



dqi 2h[ h'i 
where Q{hi) is defined by, 
'M -2-c 



hi — hi 



^ j) 



(4.25) 



Q{hi) 



aa + Bio] h^ + {{Af - 1 - c)a2 - ^n) hi + '^^ ^ ai - Bu 



hhi + bihi + 6o- 



(4.26) 



The parameters Bi^x are expressed in terms of the new three parameters bp and the five ap 

as, 

Af-2-c 



B 
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Bii = -hi + {J\f - 1 -c)a2, 
B12 



U-2 + C 
'"0 H Oi, 



^20 = 2(Ar-2 + c)a4, 
B21 = 2b2 - (SM + 1 + c)as, 
B22 = -26i + 2(Ar- 1)02, 
B02 = 2cao. 



Bqi = -26o + (A/'-2-c)ai, 
We now can easily integrate the differential equation ()4.25|) to obtain, 



(4.27a) 
(4.27b) 
(4.27c) 
(4.27d) 



>V(gi,g2) = 
where the integral constant is omitted 



1 1 + 1 1- - 1^1 - h.l (4.28) 
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E. Summary 

In summary, we have derived the canonical- form conditions for two-body case, Eqs. ()4.9p . 
and KWi . and have solved them to obtain Eqs. (HTfll . KT^ and 
Substituting Eqs. (jOD, KUl^ . KV]\i and (lOHll for Eq. (gZTD, we see that the gauged 
Hamiltonian if^ which satisfies the canonical-form condition must have the following form: 



i=l 



i=l 



_d_ 

dhi 



(4.29) 



In the above, C is calculated by substituting Eqs. ()4.18a|) and ()4.27a|) for Eq. ()3.17ep and 
reads, 



A/" — 1 

i=l 



U-l 



i=l 



where the constant R is given by, 

^_ (Ar-l)(Ar-2-c) 



2 ' hi-hj 



(4.30) 



(4.31) 



V. GL(2, K) SHAPE INVARIANCE 

It was shown that the one-body st(2) quasi-solvable models can be classified using the 
shape invariance of the Hamiltonian under the action of GL{2, K) {K = M or C) of linear 
fractional transformations 0, S|. We can see that the two-body Hamiltonian ()4.29p also 
has the same property of shape invariance. The linear fractional transformation of hi is 
introduced by, 

h^^h, = ?^l^±l {a,P,j,SeK; A = aS-Pjj^O). (5.1) 
-fhi + d 

Then, it turns out that the Hamiltonian ()4.29p is shape invariant under the following trans- 
formation (with M = 2) induced by Eq. ()5.1|) . 

M M 

H^h) ^ H^h) = Hi^h, + 6)^-' H^h) + 5)-(^-i), (5.2) 

i=l i=l 

where the polynomials P{hi) and Q{hi) in the Hj\f{h) are transformed according to, 

P{h,) ^ P{hi) = /\"\ihi + 5YP{hi), (5.3a) 
Q{hi) ^ Q{hi) = A-\^hi + dYQihi). (5.3b) 
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This shape invariance originally comes from the GL{2, K) invariance of the solvable subspace 
(ESI): 

M 



VM{h) ^ VN{h) = span | \{{ihi + 5)^-^^ : G V^ih) 



where Vj\j-{h) is given by, 



M 



V^Ch) = span <^ aiChT' ■ ■■aMihT^' ■ G Z>o, < < AT - 1 



(5.4) 



(5.5) 



From the invariance (j5.4p . the gauged Hamiltonian (|3.16|) must be invariant under the 
GL{2, K) transformation (|5.2|) because it is constructed so that Hj\fVj\f C Vv. Therefore, 
the gauged Hamiltonian ()4.29p constitutes an invariant subspace of the GL{2, K) invariant 
operators for M = 2 which satisfy the canonical-form condition. 

From the above consideration, the M-body extension (M > 2) of the two-body gauged 
Hamiltonian ()4.29|) should have the same shape invariance property under the GL{2, K) 
transformation ()5.2|) . From the functional form of the two-body sl{3) quasi-solvable system 
(j4.29|) and also from the form of the one-body sl(2) quasi-solvable system 3,, ^, ^, 22], we 
can cast the M-body quasi-solvable gauged Hamiltonian which contains up to two-body 
interactions as. 



M pa ^ a 

i=l * i=l 



dhi 



M 



^ hi — hj dhi 
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M 



C2{M)Y,P"{hi 



1=1 



N 



i=l 



hi — hj 



(5.6) 



where Ci, C2 and C3 satisfy the following initial conditions: 

U -2 U -2 + c 



Ci(l) 



C,{2) 

^2(2) =Ar-2 + 2c, 



(5.7a) 

(5.7b) 
(5.7c) 



C2(l)=Ar-2, 
C3(2) = c. 

If we make the GL{2, K) transformation on the above Hj^ according to Eq. ()5.2j) . we obtain, 
HAh) = - E ^(^*)^ + E [Ci{M)P^{h,) - Qih 



i=l 



i=l 



dhi 



M 



P{h,) d 



AT-l 



hi — hj dhi 

M 

E^'(^^) 

1=1 



Af 



12 



M 

-C2{M)Y,P2{h 



i=l 



... hi - hj 



(5.8) 
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where Pi, P2 and P3 are given by, 

Pi(/i,) =A-\^h, + 6)'P'{k) + 2 -^~^+/^^~^^^ A-%(7/i. + 5)3p(/i,), (5.9a) 

L'i{M ) 

+ 12 '^"'^^(1^''^' ^'V(7/^. + 5fP{k\ (5.9b) 

P3(/i,) = A-i(7/i, + 5fP\k) + 74^A"'7(7/i. + 5fP{k). (5.9c) 

Useful formulae needed in the above calculation are summarized in Appendix El Hence, the 
requirement that the M-body gauged Hamiltonian is shape invariant under Eqs. ()5.2j) and 
()5.3|) reads, 

Pi(/i,) = P^{hi) = P\hi), P2{hi) = P"{hi). (5.10) 
Comparing Eqs. ()5.9|) with ()A5|1 . we obtain conditions for the invariance: 

Ci(M) = ^~^^^^^~^^^ C2(M)=Ar-2 + 2(M-l)c, C3(M) = c, (5.11) 

which in turn satisfy the required initial conditions ()5.7|1 . Thus, the requirement of the 
GL{2, K) invariance uniquely determines the M-dependence of the coupling constants as 
was mentioned before. Finally, we obtain. 



M „2 r 

HM{h)=-Y,p{h;)—^+Y. 

i=l * i=l 



^-^+i^^-'>V (;.)-Q(/, 



d_ 

dhi 



/i,; — hj dh 

where C is given by, 



-2-i:Sr7^-^W"))- (5,12) 



C(a(/^)) = [AT - 2 + 2(M - l)c] 5^ P"(/^0 

i=l 

- ^ )^ g {K) - ^ ^ + R. (5.13) 

The P and Q in Eqs. ()5.12j] and ()5.13j) are again a fourth- and a second-degree polynomial, 
respectively: 

P{hi) = a^hf + a^h^ + 02^? + Oi^i + Qo, (5.14a) 
Qih,) = + hh, + bo. (5.14b) 

There are 10 parameters, namely, ap {p = 0, ... ,4), bp (p = 0, 1, 2), c, P, which characterize 
the quasi-solvable Hamiltonian (j5.12j) . The function h{q), which determine the change of 
variables, is given by a solution of the differential equation, 

h'iqY = 2P{h{q)). (5.15) 
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If we transform back the gauged Hamiltonian ()5.12j) with h{q) satisfying Eq. ()5.15|) . the 
original Hamiltonian becomes the Schrodinger type, 




dqf 



C{a{h)) 



(5.16) 



and the superpotential W(g) is given by, 



M 



Qjhj) Af-l + {M-l)c 
2P{hi) ^ 2 



M 



M 

c \n\hi — hj 

i<j 



(5.17) 



i=l 



By construction, the elements of the vector space give the solvable wave functions V'(?) 
of the Hamiltonian ()5.16|1 : 



From the form of the Hamiltonian ()5.16|) . we can observe an interesting feature of the 
(quasi-)solvability. As we have examined in the end of Section IIIIl the system is solv- 
able if Eq. 1)3.211) is fulfilled. The resulting equation ()3.23p together with the form of the 
Hamiltonian ()5.16|) indicate that the system becomes of supersymmetric form (except for 
the irrelevant constant term —C) [i^ li^ if it is not only quasi- solvable but also solvable. 
A system is always quasi-solvable if it is supersymmetric, since the ground state is always 
solvable. From the above result, we can conclude that a system is always supersymmetric 
if it is solvable and all its states have the form ()5.18j) . We note that this interesting feature 
does not depend on the normalizability of the solvable sector, that is, it holds whether or 
not supersymmetry is dynamically broken. 

VI. INVERSE PROBLEM 

In the preceding section, we have derived the M-body Hamiltonian containing up to 
the two-body interactions ()5.12|) . without directly solving the canonical- form condition for 
M > 2. To confirm that the Hamiltonian ()5.12|) is indeed quasi- solvable, i.e., it preserves 
the vector space ()3.6|) for any integer M, we shall see in this section that it can be surely 
written in terms of the elementary symmetric polynomials (Xj as the form of Eq. ()3.16|) with 
the set of the specific form of the polynomials ()3.17|) . In order to express the Hamiltonian 




(5.18) 
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()5.12p in terms of cxj, we need to know the following transformation laws from hi to (Ti. 



M 



M 



1=1 * kl=l 



dokdoi 



M 



M 



1=1 



k=l 



d 

dcTk 



(p = 0,...,4), 



(p = 0,...,3), 



M 



Y,K = B^^\a) (p = 0,l,2), 



i=l 



M 



^ hi — hi dhi 



M 



d 
dak 



(p = 0,...,4), 



h^ 



M 

^ h- - h ■ 



k=l 

&l\a) (p = 0,...,3). 



(6.1a) 
(6.1b) 
(6.1c) 
(6.1d) 
(6.1e) 



Substituting the above for Eqs. ()5.12|) and ()5.13|) . we rewrite the gauged Hamiltonian as, 



M 



k,l=l \p=0 



kl 



92 



dakdai 



M 

E 

k=l 



k 



,P=0 



Ar-2 + (M-l)c 



5^(p+l)a,+i5?) + 2cJ]a,cl^) 

p=0 p=0 



d 
dak 



C(a) 



(6.2) 



with, 



C(a) 



Af-1 
12 



[AT - 2 + 2(M - l)c] J2(P + 2)(P + l)ap+2B 

p=0 



(p) 



p=0 



c^(p+l)ap+iCi^^ (6.3) 

p=0 



In Appendix 121 we give the explicit procedure to obtain the transformation formulae ()6.ip . 

Let us begin with investigating C{a). Applying the formulae ()B9|) and ()B18|) to Eq. (j6.3|) . 
we obtain, 

C{a) =(7V- l)(Ar-2)a4(Ti2 -2(7V- 1) (TV - 2 + 0)04^2 



/A/- Ar-2-(M-l)c 
(AT - 1) 62 -a3 I f^i 



+ (^-l)^-^-(^-^)- 
6 



TV- 1 
Ma2 Mbi + R. 



(6.4) 



Comparing the result ()6.4p with Eq. ()3.17ep . we can easily see that the C{a) in the Hamil- 
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tonian ()5.12p is indeed of the form ()3.17ep with, 
Ao(a-) = a^al, 

^ A ^^-2-(M-l)c 
2^ Bioai = 2{M - 2 + c)a4(T2 + I 62 -a-^ \ (Xi, 



i=l 



2 



6 



Mbi + R. 



(6.5) 
(6.6) 

(6.7) 



These results are in agreement, in the case of M = 2, with those obtained in Section HVl 
namely, Eqs. (I4.18al) . (I4.27ajl and (IOT|l . 

In the next, we will examine the structure of the second-order derivative operator. Ap- 
plying the formulae ()B37|) to the first term of Eq. ()6.2|) . we obtain. 



p=0 

- a4(Ti(TfcO"i+i + Aki{(r), 



(6.8) 



where Aki{a) represents the part involving up to the quadratic terms in at and is given by. 



Aki{cr) = 04 



fc+2 



^(/c - m + 2)ak-m+2(yi+m + ^(A; - / - 2m + 2)ak-m+2cri+r 



m=0 



m=3 



as 



+ a2 



+ ao 



k+l 



^{k-m + l)afe_„+iai+^ + ^(/c - / - 2m + 1) 



m=0 



m=2 



fco-fcO"; + ^(A; - / - 2m)(Tfc_mO-j+»T 



m=l 



fe-1 



^(fc - / - 2m - l)(TA:-m-lO-;H 



m=0 



fc-2 



(M - Z + l)crfe_icri_i + ^(A; - / - 2m - 2)(rk-m-2(^i+m 



m=0 



(6.9) 



In Eq. ()6.9|) . we note that o"o = 1 and cr„fc = for A; > 0. Comparing Eq. ()6.8|) with Eq. ()3.16|) 
and ()3.17ap - ()3.17c|) . we can observe that the second-order derivative operator in Eq. ()5.12|) 
are expressed in the form of the first line in Eq. (|3.16p with the substitution (j6.5p and 



Afc(cT) = (204^2 - a3ai)ak + 2a4aiak+i, 
Akk{<y) = Akk{<y), 

Akiia) = Akiia) + Aikia) {k > I). 



(6.10) 
(6.11) 
(6.12) 



In the case of M = 2, these results are also in agreement with Eqs. (j4.18|) obtained in 

Section HVl 

Finally, if we apply the formulae ()B15|1 and ()B28|1 to the second term of Eq. (j6.2p . we 
can find that the first-order operator in the gauged Hamiltonian (j5.12p has the form of the 
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first term in the second line of Eq. ()3.16|) witli the following relation: 



M M 

i>j{>k) i=l 



{k + 2) [2{N' -2) + {k + l)c] 




0-3 CTfc+l 



{k-M - l){k - M - 2)caocrk-2- 



(6.13) 



From this relation, we can determine Bi^ and Bq^ uniquely since the coefficients Aijjk in the 
l.h.s. of Eq. (jOT!|l are fixed by Eqs. (jOTbj) and (jFTTTll-fjO^ . In the case of M = 2, we 
can again check that Bik and Bok calculated by Eq. ()6.13|1 are in complete agreement with 
Eqs. fl4.27bjl - (l4.27dj) derived in Section HVI 

Summarizing the analyses, we have shown that the M-body gauged Hamiltonian ()5.12|) 
can be written in terms of cxj as the form of Eq. ()3.16p . which ensures the quasi-solvability 
with the solvable sector ()3.6|) . In addition, we can find out a reason why the form of the 
quasi-solvable gauged Hamiltonian should be restricted to Eq. ()5.12|) . For example, the 
formulae ()B9|) and ()B11|) indicate. 



From the fact that -Bo(o") in Eq. ()3.16|) is a polynomial of at most second-degree, the highest 
degree term in the r.h.s. of Eq. ()6.14|) can be fit in the form of Eq. ()3.16|) only if p < 3. This 
explains why the coefficient of the first-order derivative in Eq. ()5.12|) must be a polynomial 
of at most third-degree in h^. Similar observation on the other parts of the operators can 
also explain the structure of Eq. (j5.12|) . 

Finally, we will return to the condition for the solvability. As we have examined in 
Section 11111 the operator (j3.16p becomes solvable when the condition (j3.21|) is satisfied. 
Comparing Eq. ()3.21|) with Eqs. ()6.5|) . ()6.6p and ()6.10p . we obtain the solvability condition 
in terms of the free parameters: 



Since the system has the GL{2, K) shape invariance ()5.2|) . it is solvable if the polynomials P 
and Q can be transformed according to Eqs. ()5.3|) so that the conditions ()6.15|) are satisfied. 
In the next section, we will see what kind of solvable models as well as quasi-solvable models 
can be constructed from the general form of the Hamiltonian (j5.16j) . 

VII. CLASSIFICATION OF THE MODELS 

For a given P{h), the function h{q) is determined by Eq. ()5.15|) and a particular model 
is obtained by substituting this h{q) for Eqs. ()5.13|) . ()5.16|) and ()5.17|) . The solution of 




(6.14) 



0.3 = cla = b2 = 0. 



(6.15) 
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Eq. ()5.15|) is easily obtained as an elliptic integral: 



dh 



(7.1) 



We note that the above solution has the reflection and translational symmetry; if h{q) is a 
solution of Eq. ()5.15|1 . h{—q) and h{q + qo) also satisfy Eq. ()5.15|1 . Furthermore, this elliptic 
integral also has GL{2,K) invariance discussed in Section IVl 



dh 



dh 



2p{h) 



dh 



{-ih + 6y 



2P{h) 



dh 



2P{h) 



(7.2) 



The elliptic integral ()7.ip is classifled according to the distribution of the zeros of P{h), 
e.g., multiplicity of the zeros. Since the transformation ()5.3a|) induced by Eq. ()5.1|) does not 
alter the distribution, the shape invariance ()7.2j) together with ()5.2j) enable us to classify 
the resultant quasi-solvable models (j5.16p according to the distribution of the roots of P{h). 
This idea was flrst introduced in Ref. |J| to classify the one-body s[(2) quasi- solvable models. 
Under the transformation (j5.3ap of GL(2,M) or GL(2, C), every quartic polynomial P{h) 
with real or complex coefficients is equivalent to one of the eight or five forms, respectively, 
shown in TablelH In TableHl i/, (72, gs & K according to the transformation group GL{2, K), 



Case 


GL(2,M) 


GL(2,C) 


I 


1/2 


1/2 


II 


2h 


2h 


III 


2i//i2 


2i/h^ 


iir 


i^{h^ + 1)2/2 




IV 


2i/(/i2 _ 1) 


2u{h^ - 1) 


IV' 


2i/(/i2 + 1) 




V 


2/l3 - g2h/2 - 53/2 


2h^ - g2h/2 - 53/2 


v 


i^{h^ + l)[{l-k^)h^ + l]/2 





TABLE I: The representatives of P{h) under the GL{2,M) and GL(2,C) transformations. 



and u ^ 0, < k < 1, g2 — '^^g'l 7^ 0. Since case III' under the GL(2,]R) can be obtained 
from case III by a GL{2, C) transformation, the potential in case III' is regarded as another 
real function representation of the complex potential real-valued on M*^ in case III. The 
same relation holds between case IV and IV', and between case V and V, respectively. The 
functional form of the potential is obtained by substituting the superpotential (|5.17p for 
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Eq. ()5.16|) . which reads (up to a constant term), 



V{q) 



2 ^ 

i=l 
M 



[dqj 

■AT. 



1 + (M 



C(a(/.)) 



l)c 



P'{K) - Q{h, 



X 



X 



'Af+1 + (M- l)c 



P'(/i,) - Qih, 



M 



M 



{MM)Y.hl-a,{M,M)Y,h, 



i=l 



M 



+ C C 



P{K) + P(/i, 



— 2a4hihj 



where the couphng constants 04 (M, A/") and a^{M^N) are given by, 

a4(M, AT) = [U^ + 27Vr(Af - l)c + M(M - l)c^ - (M - l)c - l] 04 
a3(M, AT) = [7V^ + 27Vr(Af - l)c + (M^ - l)c^ - 2(M 
- [Ar+(M-l)c]62. 



(7.3) 



(7.4a) 



(7.4b) 



For each case in Tabled we have specific values of the parameters Op and the transformation 
function h{q) obtained by Eq. 1)7.11) . Substituting all of them for Eq. ()7.3|) . we obtain a quasi- 
solvable potential in terms of qi in each the case. Due to the double poles at hi = hj {i 7^ j) 
in the last term of Eq. ()7.3)1 . the two-body potentials in all the cases have singularities at 
li = Qj 7^ j)- Hence, each of the models is naturally defined on a Weyl chamber. 



(0 <)qM < ■ ■ • < gi < 00, 
if the potential is non-periodic on M*^, or on a Weyl alcove, 

< gA/ < ■ ■ ■ < gi < 1], 



(7.5) 



(7.6) 

if the potential is periodic on M.^^ with a real period f2 [sOl- As we will see below, case I, II, 
III and IV*^'-* with i/ > 0, and III' with < correspond to the former while the other cases, 
namely, case III and IV^ ■* with < 0, III' with z/ > 0, and V'- •* correspond to the latter. In 
the latter case, a system can be quasi-exactly solvable unless a pole of a one-body potential 
in the system exists and is in the Weyl alcove. On the other hand, quasi-exact solvability in 
the former case depends mainly on the asymptotic behavior of Eq. ()5.18|) at |gj| ^ 00. In 
both the cases, finiteness of the contribution from qi ~ qj {i ^ j) to the norm of the wave 
functions ()5.18|) requires c > —1/2. 

In the following, we will show the functional forms of both the W(g) and V{q) without 
irrelevant constant terms. 



Case I 



In this case, oq = 1/2 and oi = 02 = 03 = 04 = 0. From Eq. ()7.1|) . we have, 

h{q) = q. (7.7) 
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The W(g) and V{q) are calculated as follows: 



, M , M M M 

>^ ( ^) = - 1 E - 1 E ^0 E - ^ E - I , ( r . 

i=l 1=1 1=1 i<j 



^ M M 



i=l i=l 

M 



This is the rational A type Inozemtsev model j51|, |5j, |53[ . Inozemtsev models are known 
as a family of deformed CS models which preserve the classical integrability. The main 
difference between quantum and classical case is that the quantum quasi-solvability holds 
only for quantized values of the parameters, namely, for integer values of M and A/", while the 
classical integrability holds for any continuous values. This is one of the common features 
quantum quasi-solvable models share. The quantum solvability is realized when 62 = 0. In 
this case, the model reduces to the rational A type CS model with the harmonic frequency 
hi- The parameter 60 corresponds to the freedom of the translational invariance and is 
irrelevant. 



B. Case II 

In this case, ai = 2 and ag = 02 = 03 = 04 = 0. From Eq. (jT.lj) . we have, 

h{q)=q\ (7.10) 
The W(g) and V{q) are calculated as follows: 

, M , M M M 

Mq) = -| E - t E - ^1 E 1- i^^i - ^ E 1- k.' - I > (7-11) 

i=l i=l i=l i<j 

^ M M 

V{q) = - E ^'(^2?.^ + h)' + -[AAf + 4(M - l)c + 2ci - 1] 62 E 



i=l i=l 
M , M 



cAci — 1) 1 , ^, ■r-^ 

+ ^T-^E^+^(^-i)E 

i=l ^* i<j 



1 1 

+ 



(7.12) 



The parameter ci is introduced according to, 

2ci = 60 -A/'+ 1 - (M- l)c. (7.13) 

This is the rational BC type Inozemtsev model. The quasi-exactly solvable model reported 
in Ref. [i^ is just this case. The solvability is realized when 62 = 0. In this case, the model 
reduces to the rational BC type CS model with the harmonic frequency 61 /2. 
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C. Case III 



In this case, 02 = 2v and = ai = = = 0. From Eq. ()7.H) . we have, 

h{q) = e^^". (7.14) 
The and V{q) are calculated as follows: 

, A/ , M M 

= - E + E - v^r^i - AT + 1) 5^ 

i=l j=l 4=1 

A/ 

- c In |sinh ^/^{qi - qj) \ , (7-15) 

i<j 



V{q) =^Y1 + i[h+Af+{M- l)c] J] e'^"^ 

j=l i=l 

, 2 Af , Af 

+ ^ V e-^^''^ + ^ [61 - AT - (M - l)c] V e-2^/^^^ 

i=l i=l 
^ 1 

This is the hyperbolic (z/ > 0) and trigonometric {u < 0) A type Inozemtsev model. The 
solvability is realized when 62 = 0. In this case, the model reduces to the hyperbolic 
(trigonometric) A type CS model in the external Morse potential. The parameter 61 = hi/2i> 
is related to the translational invariance of the model and is irrelevant. 



D. Case III' 

In this case, = i^/2, 02 = z/, 04 = z//2 and ai = 03 = 0. From Eq. ()7.1|1 . we have, 

h{q) = tan y/Dq. (7-17) 
The W(g) and V{q) are calculated as follows: 

A/ M M 

yV{q) = - 2y^C3 ^ Q'i - C2 ^ sin 207^^ - ci ^ cos 2^/l}qi 

i=l i=l 1=1 

M M 

- (A/" - 1) ^ In |cos v^gj| - c^ln \sm ^/l}{qi - qj) \ , (7-18) 

i=l i<j 
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V{q) = {4 - c?) z/^cos4v^gi + 2z/|2c2C3 - [Af + {M - l)c]ci| ^ cos2v^gi 

i=l i=l 
M M 

- 2ciC2yy^^sm.A^qi - 2z/|2ciC3 + [M + {M - l)c]c2| ^sin2yz/gi 

i=l i=l 
M 

^ sm Vi^igi - gj) 



+ z/c(c 



The parameters are introduced according to, 

hi 



Cl = - — , C2 

4z/ 



&2 - 

4i/ ' 



C3 



&2 + fep 

4z/ 



(7.20) 



The solvabihty is reahzed when C2 = C3. The parameter C3 is related to the translational 
invariance of the model and is irrelevant. As was mentioned before, this case gives another 
real function representation of case III. 



E. Case IV 

In this case, = — 2z/, 02 = 2z/ and ai = 03 = 04 = 0. From Eq. (j7.ip . we have, 

/i(g) = cosh 2v^g. (7.21) 
The W(g) and V{q) are calculated as follows: 

AI AI AI 

W(g) = — C3 cosh 2\fvqi — C2 In | sinh l^puqi | — C\ In | sinh \fvqi | 



M 



— c In I sinh \/i'{qi — qj) sinh ^/^{qi + ?j) | 



(7.22) 



M 



V(g) = 2ucl ^ sinh^ 2v^gi + 4z/C3 [ci + 2c2 + 2(M - l)c 



i=l 



M 



AI 



+ 2A/' - 1] ^ sinh^ v/^Jg^ + 2z/C2(c2 - 1) J] 



j=i 



^ sinh 2y/T'qi 



+ i/c(c-l)5^ 



+ 



sinh yz7(gi - g^) sinh ^/T'{qi + 



The parameters q are introduced according to, 

h + h2 bi-bo-b2 lr«r 1 ^2 

Cl = — , C2 = -[A/ - 1 + (M - l)cj, C3 



2i^ 



4z/ 



4i/ 



(7.23) 



(7.24) 
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This is the hyperbohc (z/ > 0) and trigonometric (z/ < 0) BC type Inozemtsev model. The 
solvabihty is reahzed when 62 = 0, that is, C3 = 0. In this case, the model reduces to the 
hyperbolic (trigonometric) BC type CS model. The quasi-solvable models constructed by 
an ansatz method in Ref. js^] are all included in case I, II, III' and IV. 



F. Case IV' 

In this case, oq = 2z/, 02 = 2z/ and ai = 03 = 04 = 0. From Eq. ()7.1|) . we have, 

/i(g) = sinh 2v^g. (7.25) 
The W(g) and V{q) are calculated as follows: 

M M M 



>V(g) = - C3 ^ sinh 2^qi - C2 ^ gd 2^/b'qi - ci ^ In [cosh 2^/uqi | 

i=l 1=1 i=l 

M 

- c ^ In I sinh v^(gi - qj) cosh y/iy{qi + qj) \ , (7.26) 



M M 

V{q) = 2ucl cosh^ 2y^qi + 2vOi [2ci + 2(M - l)c +2M -l]^ sinh 2v^g, 

i=l 2=1 

1 1 



i<j 



sinh A/z/(gi - qj) cosh v^(gi + 



(7.27) 



where gdg = arctan(sinh g) is the Gudermann function. The parameters q are introduced 
according to, 

c.-l-lW-lHM-l)c], e.^^, C3^^. (7.28) 



This form is neither the Inozemtsev nor the Olshanetsky-Perelomov type [50| even in the 
solvable case, ^2 = C3 = 0. However, it can be expressed as the one in case IV if the param- 
eters in the potential are extended to be complex. The quasi-solvable models constructed 
in Ref. [i^ are all included in case I, II, III, IV and IV'. 



G. Case V 

In this case, Oq = —(73/2, ai = —(72/2, 03 = 2 and 02 = 04 = 0. From Eq. ()7.1|) . we have, 

h{q) = p{q;g2,g3), (7.29) 
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where p{q', 92,93) is the Weierstrass eUiptic function with the invariants 92 and 93, and will 
be abbreviated to p{q) hereafter. The W(g) and V{q) are calculated as follows: 



A/ „ , X A/ „ , M 



ttJ p'iii) UJ ^ (^^) tr 

-c^ln|p(g,)-p(g,)|. (7.30) 



M M 

+ 2ci(ci - 1) 5Z P(2g^) + c(c - 1) X] - + ^(^^ + ^i)] • C^-Sl) 

i=l i<j 

The parameters 60 > Ci are given by, 

&o = &o + ci = I - ^[AT- 1 + (M - l)c], (7.32) 

and the coupling constants ai and C2(M, A/") are defined by, 

«i = ^[66o(2ci-l) + 6?], (7.33a) 

«2 = y[(2ci-l)(?2 + 26o], (7.33b) 

«3 = ^ [I86i(2ci - 1)(73 + I2hl + 52^2] , (7.33c) 

C2(M, AT) = [TV - 1 + (M - l)c + 2ci] [27V - 1 + 2(M - l)c + 2ci] . (7.34) 

The above expressions ()7.3()|) - ()7.3H1 can be arranged into a more familiar form if we introduce 
the values of p(g) at the half of the fundamental periods 2uJi: 

ei = p{uji) (/ = 1,2,3). (7.35) 

Utilizing the well-known formulae of the Weierstrass function: 

2P{h) = p'iqf = A{p{q) - d) (p(g) - 62) {p{q) - 63), (7.36a) 

2P'ih) = 2p"{q) = Upiq^ - 92, (7.36b) 



and the addition theorem, 



Jj2 

p{q^ + uji) = —-^ + ei, (7.37) 

p{qi> - ei 
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we have the following forms of yV{q) and V{q): 



3 , T M M 

«=1 ' i=l 1=1 



C 

i<j 



Y,in\p(qi) - p(qj)\ . (7.38) 



Z=l ' i=l 1=1 

M M 

+ 2ci(ci - 1) X] ^(2?i) + c(c - 1) E - + + ^i)] • C^-^^) 

i=l i<j 



M 



In the above, Hi are defined by, 
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Hf = (ei - em){ei - e„) = Se^ - ^ (/ = i, 2, 3; / ^ m ^ n ^ /). (7.40) 

The potential form (j7.39p shows that this case is the elliptic BC type Inozemtsev model. 
Contrary to all the previous cases I-IV, the solvability is not realized for any values of the 
free parameters because as 7^ 0. On the other hand, the model becomes the elliptic BC type 
CS model when ai = 0:2 = as = 0, or equivalently, 60 = ^1 = 0. Therefore, in contrast to the 
fact that the rational and hyperbolic (trigonometric) A and BC type CS models are solvable 
as quantum systems, the elliptic BC type CS model can be only quasi-solvable but is not 
solvable as a quantum system. The quantum elliptic model investigated in Ref. [s^ [H^ is 
just the CS model with the translation Qi ^ qi + We note that from the formula, 

p{qi + Wi) = pp/^^ (qi) - p{qi) + e,, (7.41) 

where p\^^'^\q) denotes the Weierstrass function with one of the fundamental periods 2uji 
replaced by ui, the above model is related to the twisted elliptic CS models 5^ 5^, 57^ . 



H. Case V 

In this case, Oq = 02 = z/(2 — k'^)/2, 04 = i/(l — k'^)/2 and ai = = 0. From 
Eq. (|7.1|) . we have, 

M«) = (7.42) 
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where sn{^/l/q\k) etc. are the Jacobian elhptic functions with the modulus k, and will be 
abbreviated to sn^y/uq) etc. hereafter. The W(g) and V{q) are calculated as follows: 



W(g) 



62 /" sn^ ^/uqi bo ^ f , cy? ^qi 



1=1 



1=1 



dn ^/lyqi 



bi Af-l + CM 

M 



V M M 

^ln|dnv^g,| - (A^- 1)^ In I cn \^q, 

^ i=l i=l 



; In |sn i/z/^i cn ^/uqj — cn ^/l'qi sn ^/l'q 



(7.43) 



i<j 



V{q) = «! > sn y/uqi cn Vz/g^ + as > sn y/uqi + } 

^-^ ^-^ ^-^ dn vz/gj 



1=1 



M 



vdc — 1) v-^ 



dn^ A/z/gj + dn^ v^g^ 



;sn ^qi cn ^/z^gj - cn ^qi sn ^/^^g 



A2' 



where cm = {M — l)c and the coupling constants are given by, 
bo - &2 



^2 



as 



04 



2A;2z/ 
A;'26o - 62 



[261 - (Ar + CM)fcV]^ 



2P 



61 - ^-j^'^' k-^u 1 ( 61 



2A;2z/ 



[26i + (Ar + CAf)A;V]^ 



2fc4i 



A/" - 1 + Cm 



(7.44) 



(7.45) 
(7.46) 

(7.47) 
, (7.48) 



where k'"^ = 1 — k^. The solvability is not realized for any values of the free parameters 
because 04 7^ 0. As was mentioned before, this case gives another real function representation 
of case V. 



VIII. POSSIBILITY BEYOND TWO-BODY INTERACTIONS 

So far, we have only concerned with up to two-body interactions. In this section, we will 
investigate what kind of many-body interactions can be constructed from our procedure. 
It can be achieved, in principle, by solving the canonical-form condition 1)4.11) for M = 
3,4,5, . . .. However, we will see that the GL{2,K) symmetry provide a powerful tool for 
this problem. 

Let us first consider the M = 3 case. From the same procedure as in Eqs. ()4.4|) - ()4.6|) . we 
obtain the transforms of the differential operators: 
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As we have seen in Section lIV] the second-order differential operators in the gauged Hamil- 
tonian must be one-body in order to satisfy the canonical-form condition. So, three-body 
operators can only exist as first- and zeroth-order differential operators. From the form of 
the gauged Hamiltonian Eq. (|3.16|) for M = 3 with Eq. (j8.ip . we can conclude that only the 
following form of the three-body operator may be added as first-order differential operators 
to the M-body gauged Hamiltonian (j5.12j) : 

.J-t^.ih,-h,)ih,-h)dhi' ^ • ^ 

where 0{hi) denotes a polynomial. On the other hand, Hj^ as a whole should be shape 
invariant under the GL{2,K) transformation ()5.2|) . From the formulae ()A2|) and ()A3|) . the 
operator ()8.2|) is transformed as, 

^ d 
+ (M - 1)(M - 2) J] A~V(7^i + ^YO{hi)— + (Oth-order ops.), (8.3) 



=1 



where (Oth-order ops.) denotes the terms which do not involve any differential operators. 
From the first term of the r.h.s. in Eq. ()8.3p . the requirement of the GL{2,K) invariance 
determines the transformation of the 0{hi): 

0{h,) ^ d{hi) = A-\^h, + 5fO{k). (8.4) 

This means that 0{h) should be a polynomial of sixth-degree. To complete the invariance, 
the second and third terms of the r.h.s. in Eq. ()8.3|) should be absorbed into the other 
two-body and one-body operators, respectively. We find that these operators should be of 
the following forms: 

since the derivatives of Oihi) transform according to Eq. ()A6jl . Similar observation tells 
us what kind of zeroth-order operators should simultaneously exist. Finally, we find the 
GL{2, K) invariant combination of the three-body operators: 

£.(3) ^ ^ 0{K) d _ M-2 ^ 0'{h,) d 

^ (h: - h .\(h. - hA rlh. 9, 2^ 



^- ~ ^i) ~ ^k) dhi 3 ^.hi- hi dhi 

30 ^ dhi Q (hi - hj){hi - hk) 

M-2^ 0"{hi) , (M-l)(M-2) .^^„,,^; 

+ m ^^^^ 



(8.6) 
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Therefore, the three-body terms should be added to the gauged Hamiltonian ()5.12p as the 
constant multiple of the combination ()8.6p . The form of the total gauged Hamiltonian reads, 



i=l 



i=l 



dhi 



(8.7) 



where Si{h) is given by, 

U -2 + {M-l)c^, 



M 



+ 93 



M 



0{K 



M 



{hi - hj){hi - hk) 



M-2 ^ O'jh,) (M-l)(M-2) ^„^^ 
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(8.8) 



In order that the total gauged Hamiltonian ()8.7p can be cast in the Schrodinger form by a 
gauge transformation, we must have. 



dW _ h'l S,{h) 
dqi 2h'. h[ 



M), 



(8.9) 



in addition to Eq. ()5.15|) . Thus, the integrability condition for W(q') becomes. 



d dW d dW 



h'ids,{h) h[dSi{h) 



h\ dhi 



(Vz^/). 



dqi dqi dqi dqi h[ dhi 

However, this condition cannot be satisfied unless (73 = since we obtain from Eq. 



(8.10) 



h[ dS,{h) 
K dh. 



h'M 



_93 K 

{K-hif • iK-hi^h'^ 



+ 



M 



20{h,) 
hi — hi 



M-2 



0'{h,) 



(8.11) 



In other words, the 5[(M+1) quasi-solvable operator ()3.16|) which satisfies the canonical-form 
condition cannot contain the three-body operators. 

The generalization of the above argument to the M > 3 cases is straightforward and re- 
sults in the same conclusion that the existence of M-body operators for M > 3 is prohibited. 
Therefore, the quasi-solvable quantum many-body systems constructed from the procedure 
in Section IIIII are exhausted by the ones classified in the preceding section. 



IX. CONCLUDING REMARKS 

In conclusion, we have proposed a systematic method to construct quasi-solvable quantum 
many-body systems. In Table ITTl we summarize the complete list of the quasi-solvable models 
as well as their special subclass of the solvable models constructed in this paper. 

We would like to close this paper by giving several remarks on the future problems. 
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Uase 


Quasi-solvable 


boivabie 


T 
i 


rational A Inozemtsev 


rational A Calogero-Sutherland 


II 


rational BC Inozemtsev 


rational BC Calogero-Sutherland 


TTT 
ill 


hyp. (trig.) A Inozemtsev 


hyp. (trig.) A Calogero-Sutherland 
+ external Morse potential 


IV 


hyp. (trig.) BC Inozemtsev 


hyp. (trig.) BC Calogero-Sutherland 


V 


elliptic BC Inozemtsev 
(elliptic BC Calogero-Sutherland) 





TABLE II: Classification of the (quasi-)solvable quantum many-body systems. 



1. We should stress that quasi-solvable quantum many-body models which can be con- 
structed from s[(M -|- 1) generators are not limited to the ones presented here. This 
is because it depends on the way of changing of the variables and on the choice of the 
solvable subspace Vj\f. In Ref. ji^l, several M-body quasi-solvable quantum models 
with up to M- and up to 6-body interactions were constructed from s[(M + 1) gen- 
erators. The Calogero-Marchioro-Wolfes model |^ or the rational G2 type CS 
model and the trigonometric G2 type CS model without the two-body interaction, all 
of which have three-body interactions, were shown to have sl(3) algebraizations and 
thus to be solvable in Refs. [62, 63, 64\. Therefore, there may be another family of 
quasi-solvable quantum many-body systems which can be constructed from a different 
s[(M + 1) algebraization. 

2. From Table ITTl we notice that, among the classical type CS models, only the elliptic A 
type CS model cannot be obtained from the present procedure. As far as we know, the 
quantum (quasi-)solvability of this model has not been confirmed yet. We note that 
this model is also a special in the sense that any deformation of the model preserving 
the integrability have not been discovered yet in contrast to the other classical type 
CS models. 

3. Table |n] indicates that classical integrability may have more intimate relation with 
quantum quasi-solvability rather than with quantum solvability and quasi-exact solv- 
ability. The latter is understood because classical integrability does not depend on 
whether the particles are moving inside a bounded region or not. So, it will be an 
interesting problem to study the relation between classical integrability and quantum 
quasi-solvability. 

4. Although the paper only deals with scalar problems without any internal degree of 
freedom, there have been much progress in constructing quasi-solvable spin systems 
recently. Among them, the generalized Dunkl operator approaches have provided one 
of the most powerful tools for the issues |6^ |6^ |6^. Then, it is quite interesting 
to note that the spin CS models constructed from the AM-type Dunkl operators [65| 
coincide, in the scalar limit, with the (scalar) models constructed from the s[(M + 1) 
generators in this paper. Thus, there should be some intimate relations underlying 
between the two approaches although the algebraic structure of the y^M-type Dunkl 
operators and that of the sl(M + 1) are quite different from each other. Another 
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interesting problem is the possibility to generalize the Dunkl operators to include M- 
body (M > 3) operators. The conclusion in Section IVIIII however indicates little 
possibility for such a generalization, at least, of the ^M-type. 

5. In this paper, we imposed the canonical-form condition (j4.1|) in order to obtain the 
ordinary Schrodinger operators. However, it may be possible to obtain another class of 
quasi-solvable second-order operators by imposing another type of condition instead 
of Eq. ()4.H) . If we can solve the canonical- form condition on a specific non-trivial 
metric or on a particular manifold, we will obtain a gravitational deformation of the 
quasi-solvable quantum models. 

6. A generalization of the Bender-Dunne polynomials to the models obtained here will 
be an interesting problem. Although It was argued in Ref. [3] that they exist for any 
quasi-solvable quantum systems regardless of the number of particles, we have not 
appreciated whether their arguments can be naively applicable to our case. 

7. As we have discussed in the end of Section El all the solvable models constructed in 
this paper are supersymmetric. It was shown in Refs. [i^ IH, 6^ that the rational 



and hyperbolic (trigonometric) A and BC type CS models are of the supersymmetric 
form. Table IH] shows that the hyperbolic (trigonometric) A type CS model can be 
deformed without destroying both the solvability and supersymmetry by adding the 
external Morse potential. Then, a natural question arises whether there is a solvable 
quantum model which is not of the supersymmetric form. 

In the case of M = 1, it was shown that quasi-solvability and A/"- fold supersymmetry 
is essentially equivalent However, the relation between them has not been under- 
stood in many-body quantum systems although A/'-fold supersymmetry in many-body 
case was briefly formulated in Ref. [16|]. We expect that they also have the similar 
intimate relation to each other in the many-body case. 
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APPENDIX A: FORMULAE FOR GL{2, K) TRANSFORMATION 



Here we will list the formulae useful for the calculation of GL{2, K) transformation in 
Section |Vl Under the linear fractional transformation ()5.1|) . we have. 



hi 



hi — hi 



hj 



hj 



hi — h 



' {'yh, + 5)i^hj + 5) 
1 



A-'{^h, + 6) 



hi — hj 



^{^hi + 6) 



-1 



(Al) 
(A2) 
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The differential operators in the gauged Hamiltonian are transformed according to, 

M M 



_d_ 

dhi 



iN-l) 



_d_ 

dhi 



(AT- 1)7(7/1, + 



(A3) 



and, 



M 



M 



J.V1 rj2 
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-(A^-l) 



+ {U-l){N-2)^\^hi + 5)- 



(A4) 



From the transformation rule (j5.3|) for the polynomials P{h) and Q{h), the derivatives of 
the transformed polynomials are calculated as, 

P'{h,) = A-\-fh, + 5fP'{U) + 4A-27(7/i, + 5fP{U), (A5a) 
P"{hi) = P"{k) + 6A-S(7/^. + 5)P'{k) + 12A-V(7/^. + 5fP(hi), (A5b) 
Q\K) = Q'{k) + 2A-S(7/i, + 5)g(/i,). (A5c) 

Similarly, the derivatives of the transformed sextic polynomial 0{h) given by Eq. ()8.4|) are 
calculated as, 



d\hi) = A-2(7/i, + 5fO\k) + 6A-^7(7/i, + <5)^0(/i,), 
6"(/i,) = A-i(7/i, + 5)'0"(/i,) + 10A"-%(7/i, + 5)30'(/i,) 

+ 30A-V(7/^^ + ^)'O(/i,), 
6'"(/i,) = 0"\k) + 12A-S(7/i. + 5)0" [k-:) 

+ 60A-V(7/ii + ^)'0'(/ii) + 120A-Y(7/^i + 5)'0(/ii 



(A6a) 
(A6b) 
(A6c) 



APPENDIX B: TRANSFORMATION FORMULAE BETWEEN hi AND ct. 

In this appendix, we will show a systematic method to calculate the transformation 
formulae between the variables hi and cxj needed in Section IVII and IVIIII The method 
presented here is essentially based on the one in Ref. ji^l. The starting point is that, for the 
elementary symmetric polynomials cxfc and the Newton polynomials Sn defined by, 

M 

(^k{h) = £^0 = 1; s„(/i) = ^/i^, (Bl) 
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the following relation holds: 

G{t) = J2 = exp 5^ M^"""^" ■ ^^^^ 

k=0 \n=l ) 

Thus, Gif) gives the generating function for both the elementary symmetric polynomials 
and the Newton ones. The derivatives of G(t) with respect to hi and t are, 

^G{t) = G{t)Y}.-r^'K-'t- , (B3a) 

* n=l 

^G{t) = G{t) = J2 k^kt'-\ (B3b) 

n=l fc=0 
pn / oo CO \ 

" G{t) = G{t){ Yl (-r^™'^'t"''™'"'sms™'-$^(-r+'^""'"^w) • (B3c) 

\m,m'=l m=l / 



dt 



We will see that most of the formulae involving differential operators can be derived with the 
aid of the generating function G{t). Here we derive a useful formula for the later purposes. 
Let us define a generating function B(t) as, 

k=0 * 

From Eq. ()B3aj) . we have another expression for B{t): 

B{t) = ^G{t) = tG{t) - Y,{-rhr't-G{t) 

* n=2 

Comparing Eq. ()B4|) and the last expression in Eq. ()B5|) . we have the following formula, 
which facilitates the derivation of the several formulae. 



1. Calculation of ^[^^(cr) 



The definition of the quantity B^^ is given by, 

M 



J2K = s,{h)=B^'\a). (B7) 



i=l 
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From the formula ()B3b|) . we can immediately obtain 



The first few formulae read, 



t=o 



(B8) 



B 
B 
B 
B 



(0) 

(1) 

(2) 



(3) 



so = M, 



Si = CTi, 



52 = O"! - 2(72, 

53 = O"! - 3(72Cri + 30-3. 



(B9a) 
(B9b) 
(B9c) 
(B9d) 



2. Calculation of sl^Vcr 



The definition of the quantity B^^^ is given by, 



M 



j=l A:=l 



B 



ip) 



(BIO) 



The usefulness of the formula ()B6|) is understood if we note that the following recursion 
relation can be derived from it: 



M 



- B_ Ok- 



The above relation enable us to calculate all the B'f^ needed, once one of them is known. 
The calculation of B^^^ will be the easiest. Let us define a generating function B^^\t) as. 



(Bll) 



M 



B^^\t) = Y^Bih^ 



k=0 



From the formulae ()B3|) . we have another expression for B^^\t): 

M r% 00 

o 



i=l 



n=l 



d 



M 



t-G{t) =Y,kakt\ 



k=0 



By comparing Eq. ()B12|) with the last expression in Eq. ()B13|) . we obtain. 



B 



(1) 



kcTk, 



(B12) 



(B13) 



(B14) 
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Applying the recursion relation ()B11|) to Eq. ()B14|) successively, we finally get the formulae 



for b'^^: 



= {M-k + l)afe_i, (B15a) 

b\^^ = kak, (B15b) 

= c^ifXfc-(A; + l)afc+i, (B15c) 

= {al - 2(T2)o-fe - aiCTfc+i + {k + 2)(Tfc+2. (B15d) 

3. Calculation of ci^V) 
The quantity C^^ is defined by, 

^ hP 

A modification of Eq. ()B16jl reads, 

^ M , p _ , p M p-1 

ly^j i^j 1 = 



.1=0 

With the aid of the formulae ()B9|) we easily obtain, for instance, 

= 0, (B18a) 

Ci') = i^(M-l), (B18b) 

Ci^) = (M - (B18c) 
Ci^) = (M - l)al - (2M - 3)a2. (BlSd) 

4. Calculation of C^^V) 
The quantity C^^^ is defined by, 

\^ ^ = V Mp) ^ ■ r^(P) ^ V ^» ^^fc mi ON 

^hi-h.dhi ^ dak' ^hi-h.dK ^ ' 

i^j ■' k=l i^j ■' 

In order to calculate C^'^ , we define a generating function C^^^ [t] as, 

M 

C^\t) = y^&^hK (B20) 

fc=0 
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From the formula ()B3a|l we have, 



M oo p+m— 2 

m=l 1=0 



(B21) 



On the other hand, we obtain for g > from the formula ()B3c|) . 



oo m—q+1 



m+l^m+g— 1 



mSm+1 



m=q 1=1 

oo 



m=l 



m—q+2 
m=q /=0 



- G{t) ^(-)"+H'^+'?-i(2M + m)sm+i. 

m=l 

The coefficient of G{t) in the second line of Eq. ()B22|) is further deformed as, 

oo m—q+2 M oo m—q+2 

J2^_yn-q+2^m ^ ...^^^^-Z = ( J] h^hj-^^'-' 

m=q 1=0 i^j m=l 1=0 

oo 

+ E(-)'"+H™+^-i(m + 2)sm+i + f^{t), 

m=l 



(B22) 



(B23) 



where fq(t) is given by. 



{2M-3)s2 + sl q = 0, 

q—l m—q+2 

J2^_-^rn-q+2^m (s„_,+2 ' Q > 1- 

. m=l 1=0 



(B24) 



Substituting Eq. ()B23|) for Eq. ()B22|) and applying the formula ()B3b|l . we obtain. 



M oo 



m—q+2 



m—q+2— I 



G{t) 



i^j m=l 



1=0 



- 2(M - ir^' (-^ + G{t) + /f (t)G(t). (B25) 



Combining Eqs. ()B20jl and ()B25|1 . we get another expression for G^^^'it) (0 < p < 4): 



2G^P\t) 



G(() 



(B26) 



k=0 
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Comparison of Eq. ()B20|) with Eq. ()B26|) enable us to calculate C^^ for < p < 4. The 
recursion relation derived from Eq. ()B6|) in this case is, 

= afeCi^) - Ci'l,. (B27) 
The calculation of Cj^'^ may be the easiest. From Eqs. ()B26|) and ()B27|) . we obtain the 



formulae for C, 



(p). 



k 



2Cf ^ = - (k - M - l){k - M - 2)ak-2, (B28a) 

2Cj^^ ={k- M){k- M -l)ak-u (B28b) 

2Cf ^ = -k{k-2M + l)(Tfc, (B28c) 

2Cf ^ =2(M - l)aiorfc + {k + l){k-2M + 2)ak+i, (B28d) 

2C^^^ = 2(M - l)aiVfc - 2(2M - 3)a2(Tfc - 2(Af - l)aiak+i 

-{k + 2){k- 2M + 3)(Tfc+2. (B28e) 

5. Calculation of Ag(o-) 

The definition of the quantity A^^j is given by, 

M r.2 M „2 M r. 

h " ' ^^^^^^ ' ' " ■ 

In order to calculate A^^], we introduce a generating function y4(P^(t,u) defined by, 

M 

A(^)(t,«)= J]AgtV. (B30) 

fc,Z=0 

Then, it immediately reads, 

oo 

= E (-)"^'+"'+2s^'+„.+p-2r'M"'G(t)G'(M). (B31) 



m'n'=l 



If we make the substitution m = m' + n' and n = m! — n' , the expression of y4*-^^(t,M) is 
rewritten as, 

oo m-2 / y_ \ n/2 

A^^\t,u) = Y,{-rs^^,.,{tur" (-) ^w^h 

m=2 n=2— m ^ 

oo 



t — u 

m=l 
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On the other hand, utihzing the formula ()B3b|) we have for p > 0, 



m=l 



(B33) 



where fp(t) is given by, 



m=l 



Mt, 

p-i 

1 m=l 



p = 0, 



\m+l±m—p+l 



t'^-^-^'s^, p>l. 



(B34) 



From Eqs. (jB32|l and (jB33jl . we obtain another expression for A^'''\t, 



u : 



u 



\p+i. 



tu 
t — u 



t^-P— _ u^-P— - f^(t) + f^(u) 

dt du + 



M 



t — U 



G{t)G{u) 



(B35) 



fcJ=0 



Comparison of Eq. ()B30|) with Eq. ()B35|) enable us to calculate A^l^\. The recursion relation 
derived from the formula ()B6|) in this case is. 



Atr^ = Bt'a, - Atl, = a,B[^^ - Ag,, (B36) 
The calculation of may be the easiest. From Eqs. ()B35|) and ()B36|) we obtain. 



k~2 



A 



(o: 



A 



(1 



A 



(2: 



^(3: 

A 



(B37a) 
(B37b) 
(B37c) 



■■{M -1 + l)ok-iOi-i + ^(/c - / - 2m - 2)ak-m-2cri+m, 

k-l 

■ - ^ (A; - / - 2m - l)crk-m-i(^i+m, 

m=0 

k 

■ kcTkOi + Y{k-l- 2m)ak-rn(^l+m, 

m=l 

1 k+1 

■ cTiCTfcO-/ - ^ (A; - m + l)ak-m+i(^i+m - ^ (A: - / - 2m + l)ak-m+i(^i+m, (B37d) 

m=0 m=2 

^ {(tI - 2^2) ak(Jl — CFi {(Tk+l(Tl + <Jk(Tl+l) 
2 k+2 

+ Y{k-m + 2)ak^m+2(yi+m + ^(A; - Z - 2m + 2) (B37e) 



m=0 



m=3 
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